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Abstract 

The ATLAS collaboration recently reported a 3cr excess in the leptonic -Z -\-jets + E™ lss channel. 
We intend to interpret this excess by squark pair production in the Next-to-Minimal Supersym¬ 
metric Standard Model (NMSSM). The decay chain we employ is q —* q \2 —> Q'xfZ, where x 1 
and X 2 denote the lightest and the next-to-lightest neutralinos with singlino and bino as their 
dominant components respectively. Our simulations indicate that after considering the constraints 
from the ATLAS searches for jets + E'™ ss signal the central value of the excess can be obtained 
for mq < 1.2TeV, and if the constraint from the CMS on -Z search is further considered, more than 
10 signal events are still attainable for niq < 750GeV. Compared with the interpretation by gluino 
pair production, the squark explanation allows for a significantly wider range of mq as well as a 
less compressed SUSY mass spectrum. We also show that the squark explanation will be readily 
tested at the initial stage of the 14 TeV LHC. 

PACS numbers: 
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I. INTRODUCTION 


Since the discovery of the Higgs-like particle by the ATLAS and CMS collaborations at 
the Large Hadron Collider (LHC) in 2012 [T], the main task of the LHC program has shifted 
to the searches for new physics beyond the Standard Model (SM). These searches cover a 
wide range of possible signatures, notably various combinations of jets (with or without b- 
tagging), the missing transverse energy (A™ ss ) and/or leptons. In this direction, the ATLAS 
collaboration has recently reported an intriguing excess at 3a significance in the leptonic- 
Z + jets + E™ lss channel^. Based on the full dataset in 2012, which corresponds to about 
20.3 fb^ 1 integrated luminosity at the 8 TeV LHC, the collaboration observed 29 events for 
the on-Z electron and muon pair channels in contrast with the expected SM background 
number 10.6 ± 3.2, and no excess over the SM background was observed in any other signal 
region (SR)0. 

So far there have been several attempts to explain the excess by the production of new 
physics particles, which are needed to decay with a sizeable rate into jets as well as at least 
one Z boson and one invisible particle p 112] . In the context of supersymmetric theories 
(SUSY), in order to provide enough events after the rather tight cuts employed in [2j, the 
production of a pair of moderately light gluinos was usually utilized for the explanation^ 13- 
El HI 12]. The key point in doing this is to choose a lightest SUSY particle (LSP) with 
relatively suppressed couplings to squarks, so that the gluino prefers to decay first into the 
neutralino other than the LSP, and the neutralino subsequently decays into the LSP plus a 
Z boson. In the Minimal Supersymmetric standard model (MSSM), a higgsino-dominated 
neutralino has very weak couplings to light flavor squarks, so one can naturally imagine a 
higgsino-dominated LSP and assume the decay chain g — >■ qqx ° -A qqZj(i (q represents a 
light flavor quark, and denotes a gaugino-dominated neutralino) to interpret the excess 1 . 
For this case, we note that the measured dark matter relic density is not easy to obtain if 
only the neutralino serves as the dark matter candidate (see for example, Fig.l in [T3j). In 
the Next-to Minimal Supersymmetric Standard Model (NMSSM)pUQ|, however, a singlino- 
dominated neutralino also has the property, and meanwhile if it acts as the LSP, the correct 
relic density can be achieved by multiple annihilation channels[TU [15j- So in this work we 

1 In this case, the higgsino-dominated LSP plays the same role as the gravitino in the ATLAS report [2] to 
interpret the excess. 
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are interested in the interpretation of the Z-excess in the NMSSM with a singlino-dominated 
LSP. 

In the framework of the NMSSM with a singlino-dominated LSP, the gluino pair pro¬ 
duction with the three-body decay g —* qq \2 —>• qqZ\ ? has been studied for the Z-excess 
in 0E]. These works indicated that only with simple assumptions on the relevant model 
parameters can the NMSSM explain the excess quite well. Explicitly speaking, it was found 
that after considering the constraints from the ATLAS searches for jets + E™ lss signal, the 
NMSSM can reproduce the central value of the excess, and even if one further considers the 
constraint from the CMS search for the leptonic-Z + jets + E™ lss channel which observed 
no excess in all SRs, the event number of the ATLAS on -Z signal can still reach 11, which 
is about 1.2a away from the measured central value p]. Moreover, as illustrated in [5] the 
gluino explanation can reproduce well various distributions of the excess presented by the 
ATLAS collaboration. Albeit these advantages, we still think that it is necessary to seek for 
other explanations since in the gluino explanation, the gluino mass is limited in a narrow 
range, and at same time the relevant sparticle mass spectrum must be rather compressed to 
escape the constraints (see Fig.2 of |[8j). In this work, we consider squark pair production 
as an explanation of the excess, and in order to compare it with the gluino explanation, 
we make similar assumptions on the model parameters to those of [8]. We find that in the 
squark explanation the range of the squark mass can be significantly extended in comparison 
with the gluino explanation, and moreover the relevant sparticle mass spectrum may become 
less compressed. We also find that, just like the gluino explanation, the distributions of the 
excess can also be reproduced well in the squark explanation. 

This work is organized as follows. In Section [TT] we briefly introduce our scenario for the 


excess. In Section III, we perform a comprehensive analysis of the relevant parameter space 


and present the results of our simulations on the Z-peaked excess. In Section [TV] we choose 
some representative parameter points, and exhibit their predictions on various distributions 
of the excess in comparison with the corresponding data provided by the ATLAS collabo¬ 
ration. In section [VJ we briefly discuss future test of our scenario at the 14 TeV LHC. At 


last, we draw our conclusions in Section VI 
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II. OUR SCENARIO FOR THE Z-EXCESS 


As one of the most economical extensions of the MSSM, the NMSSM contains one gauge 
singlet Higgs superfield S in its matter content. The superpotential of the general NMSSM 
is given b Y mm 

Wnmssm = Wmssm + A SH U ■ Hd + C,fS + + g S 3 , ( 1 ) 


where Wmssm is the superpotential of the MSSM without the /i term, H u and Hd are the 
SU(2)l doublet superfields, k and A are dimensionless coefficients, parameterizes the 
tadpole term and //' is a supersymmetric mass. 

In this framework, the fermionic component field of S which is usually called singlino 
S will mix with the gauginos and the higgsinos of the MSSM to form neutralinos. In the 
basis (0i = —f5,0 2 = — iW°,ijj 3 = H®, 04 = 5°,0 5 = S ), the corresponding mass matrix 
is given by m 
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where M\ and M 2 are soft gaugino masses, v u = r?sin/3 and Vd = vcos/3 are vacuum 
expectation values (vev) of the Higgs fields H u and Hd respectively, /j, e ff = /i + As with 
s denoting the vev of the singlet scalar field S, and c w = cos 6^- This matrix can be 
diagonalized by a 5 x 5 unitary matrix N, and consequently neutralinos as mass eigenstates 
are defined by 
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X° = 

3 = 1 

where the mass order m^o < m^o < ■ ■ ■ < m^o is assumed. Obviously, the matrix element 
Nij measures the size of the 0j component in state, and for the singlino-dominated 
and bino-dominated neutralinos, their masses are mainly determined by the combination 
2 ks + // and Mi respectively. Moreover, with the help of one can get the interactions of 
the neutralinos. As shown in |8], the qx®q coupling with q denoting a light flavor quark is 
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FIG. 1: The cross section of the squark pair production at the 8 TeV LHC for different gluino 
mass. 

determined by the gaugino components of ,\f, and the XiXqZ coupling is determined by the 
higgsino components of the neutralinos. By contrast, the x^Xjh coupling with h denoting 
the SM-like Higgs boson depends on all components of the neutralinos, and there may exist 
cancelations among different contributions. These characters are helpful to understand our 
explanation of the Z-excess. 

In the following, we intend to interpret the ATLAS on-Z excess by squark pair production. 
To make our explanation as simple as possible, we have the following assumptions: 

• Only the first and second generation squarks are responsible for the excess. In our 
analysis, we assume a common mass rriq for the squarks, then the cross section for the 
squark pair production depends only on rriq and rrig. We calculate the cross section at 
the NLO with the code Prospino m, and show its dependence on rriq at the 8 TeV 
LHC in Figjl] 

• The lcptonic -Z+jets+E™ lss signal is generated by the cascade decay q —> qx ? 2 ~► QXi^- 
In order to maximize this signal rate, we require both Br(q —> qx °) and Br(x 2 —> X°\Z) 
at roughly 100%, where the former requirement can be satisfied if Xi an d X2 are 
singlino-dominated and bino-dominated respectively, and only the two particles in the 
neutralino and chargino sector are lighter than the squarks, and the latter condition 
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can be realized if rriz < m^o — m^o < rrih or if the XiXjh interaction is significantly 
suppressed (see above discussion and also our previous work 0 )- 

• With the above assumptions, the parameters involved in our explanation are rrig, rrig , 

Ami = — m xl an d = m x 2 ~~ m xi ' ^ ri our discussion, we vary these parameters 

freely, but noting that the process pp —> qq —>■ XiZqx\Zq can also generate multi-jets 
_\_Ejp iss signal, we limit these parameters by the ATLAS searches for the multi-jets 
signal, which were presented in [18, [T9]. We also consider the CMS search for the 
leptonic-Z + jets + E™ lss signal [2CJj as an alternative constraint on the parameters. 

About our scenario for the Z-excess, we have following additional remarks: 

• We ad hoc require that only Xi an d X® am ong the neutralinos are lighter than the 
squarks. This will simplify our analysis, but on the other hand, since the rate and 
various kinematic distributions of the process pp —>■ qq —>■ XiZqX\Zq are decided by 
few parameters, the capability of our scenario to interpret the excess is limited given 
that the scenario must satisfy the constraints mentioned above. In fact, as implied by 
the results of [f0|, allowing the squarks to decay in multiple ways facilitates SUSY to 
balance the ATLAS signal and the constraints, and is thus able to explain the excess 
in a better way. This, however, needs an intensive scan over a higher dimensional 
SUSY parameter space, and for each parameter point, simulations on various SUSY 
signals have to be done to compare with the corresponding experimental data. Such 
calculations are very time-consuming, and are beyond the capability of our cluster. 

• Again for the sake of simplicity, we do not consider the effect of the third generation 
squarks in our analysis. These squarks have been tightly limited by the SM-likc Higgs 
boson mass, and are preferred to be heavy [2TH23] . For some optimized points in 
Figj2]for the Z-excess, we once included their contributions to the leptonic-Z +jets + 
E™ lss and the multi-jets +E™ iSS signals by assuming the degeneracy of all squarks. 
However, we did not find any improvement on our explanation due to the constraints 
we considered. 

• The assumptions on the properties of the LSP and NLSP in this work are same as 
those of our previous work [5j, where gluino pair production was used for the excess. 
This enables us to compare directly the two explanations. 
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III. Z-PEAKED EXCESS IN OUR SCENARIO 
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FIG. 2: Constant contours of the event number for the ATLAS leptonic-Z +jets + M^ iss signal 
on the Ami — A m 2 plane in heavy gluino case. For each rriq, the region between the contour 
marked by the number 12.1 and that by 24.7 can explain the excess at lcr level, and that between 
the lines marked by 5.8 and 31 respectively can account for the excess at 2cr level. The parameter 
spaces that coincide with different SUSY searches are also presented, which are right bounded 
by different types of lines. The dotted line and solid line are the boundaries coming from the 
ATLAS preliminary and updated searches for 2 ~ 6 jets + E™ lss signal respectively, and the dash 
dotted line is for the CMS constraint. For the case rriq = 700GeU, the constraint from the ATLAS 
preliminary search for 2 ~ 6 jets + Etp iss signal is too weak to be drawn on the plane, and for the 
case m,q = 970GeU, there are actually no boundaries on the plane. 


From the ATLAS analysis on the leptonic-Z + jets + E™ lss channel presented in [5], one 
can infer that the event number of the excess is 18.4 ± 6.3 after including the statistical 
and systematic uncertainties [8j. This means that, if one wants to explain the excess at 
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lcr and 2 a levels, the SUSY signal number after cuts should satisfy 12.1 < Nu < 24.7 and 
5.8 < Nu < 31 respectively. In order to hnd the parameter space that can produce the 
required event number, we fix nig = 4.5TeV (heavy gluino case) and rrig = 1.5TeV (light 
gluino case) separately, and for each case, we perform a grid scan over the parameters Ami 
and Am 2 by choosing a series of rriq. For each parameter point, we calculate the squark pair 
production rate at the 8 TeV LHC by the package Prospino ra. and generate the parton 
level events for the considered process with MG5_aMC|24j. which includes Pythia[25j for 
parton showering and hadronization. Then we use the package CheckMATE-1.2.0 [2B] which 
contains kne-tuned fast detector simulation code Delphes3.0.10[27] to repeat the analyses of 
various experiments. These experiments include the ATLAS on-Z search[2], the CMS on-Z 
search[20], and the ATLAS 2 ~ 6 jet -\-E™ lss searches 1.8. 19], among which the first one is 
used to generate the signal of the excess, and the other ones sever as constraints. In [8], we 
encoded the cuts for those experiments in the package CheckMATE-1.2.0, and the validation 
of them indicated that our calculations coincide with the corresponding experimental anal¬ 
yses at 20% level. In implementing the constraints of the SUSY searches on the parameters, 
we define for each search the ratio R = max(Ns,i/S^^), where Nsj is the event number of 
the SUSY signal in the All SR of the search, is its 95% upper limit usually provided 
in the experimental report, and the max is over all SRs defined in the search. Obviously, 
only in case that R < 1, the corresponding parameter point is experimentally allowed at 
95% C.L.. 

In Figj2]and Figj3j we present on Ami — Am 2 plane the constant contours of the event 
number for the ATLAS on-Z analysis in the heavy gluino and light gluino cases respectively. 
For each m,, the region between the contour marked by the number 12.1 and that by 
24.7 can explain the excess at la level, and that between the contours marked by 5.8 and 
31 respectively can account for the excess at 2 a level. The parameter spaces that coincide 
with various SUSY searches are also presented, which are right bounded by different types of 
lines (note that a compressed SUSY mass spectrum is helpful to evade the LHC constraints). 
The dotted line and solid line are the boundaries coming from the ATLAS preliminary and 
updated searches for 2 ~ 6 jets + E™ lss signal respectively, and they are obtained by setting 
the corresponding R values at 1. The constraint from the CMS on -Z search is obtained in 
a similar way, and shown by dash dotted lines. 

From Fig]2]for the heavy gluino case, one can learn following facts: 
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FIG. 3: Similar to Figj2j but for the light gluino case, rrig = 1.5TeV. Note that for rriq = 1360GeV, 
the CMS constraint is too weak to be drawn on the plane. 

• For all the choices of rriq, the strongest constraint comes from the CMS dedicated on-Z 
counting experiment, and the weakest one is the ATLAS preliminary search for 2 ~ 6 
jets + E™ lss signal. 

• With the increase of the squark mass, Ami is allowed to vary within a wider range. In 
this case, the improved cut efficiency due to the enlarged Ami can compensate for the 
decrease of the squark pair production rate. As a result, even for rriq < 800GeV the 
central value of the excess (18.4 events) can still be obtained if only the constraints 
from the ATLAS searches for the jets + E™ lss signal are considered. 

• In the CMS dedicated on-Z counting experiment, six signal regions discriminated by 
the jet number rij and E 7 ™ ss were considered (see Table 2 in [8]). We checked that 
for our scenario the tightest constraint of the experiment comes from nj > 2 SRs 
with Etp iss either satisfying 200GeV < E™ lss < 300GeV (called SR-II hereafter) or 
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TABLE I: Details about the constraints in light gluino case for two points Si and S 2 , which are 
taken from the lower left and right panels in Figj3] respectively, and are defined by Si: niq = 
HOOGeV, Ami = 458GeV, Am 2 = 105GeV and S 2 : Trig = 1360GeV,Ami = 687GeV, Am 2 = 
105GeV. In this table, SR max represents the SR which provides the strongest constraint on the 
parameter for a certain experiment, and e is the total cut efficiency of the SR. 


Point 

cr(fb) 

CMS HU 

ATLAS (updated) [18J 

ATLAS (preliminary) [19j 

SRmatc 

e 

c95% 

D obs 

R 

SR/natc 

e 

c95% 
^obs 

R 

SR ■max 

e 

c95% 
^ obs 

R 

Si 

32.9 

SR-III 

3.3% 

7.6 

2.77 

2jm 

14.7% 

90 

1.09 

CM 

9.2% 

81.2 

0.75 

s 2 

6.3 

SR-III 

4.2% 

7.6 

0.67 

2jt 

21.0% 

38 

0.71 

BM 

12.4% 

14.9 

1.06 


satisfying E™ lss > 300GeV (called SR-III hereafter) 2 . In either case, the signal of the 
ATLAS on -Z search has a large overlap with that of the CMS on -Z search, so due 
to the tension of the two search results the event number of the ATLAS experiment 
is always upper bounded by about 11 for rriq < 750GeV after considering the CMS 
constraint. We also checked that, with the further increase of rriq from about 750GeV, 
the maximal reach of the event number drops either because the tension of the ATLAS 
and CMS data becomes stronger for moderately heavy q or because the squark pair 
production rate is sufficiently suppressed for heavy q. 

• The lower right panel of Figj2] indicates that there are actually no boundaries on the 
Am! — Am 2 plane for rriq = 970GeV. In this case, the maximal reach of the event 
number is 6.8, which is still within the 2 a range of the excess. 

Next we turn to the light gluino case. From the event contours in Fig|3j one can get 
following information: 

• The dependence of the ATLAS event number on the parameters A mi and A m 2 is 
quite similar to that in the heavy gluino case, and so are the dependencies of the 
constraints. 

2 In more detail, from our calculations we learn that the SR-II is more powerful than the SR-III in limiting 
our scenario for the region defined by lOOGeV < Ami 1$ 150GeV and 95GeV < Am 2 < 150GeV as 
well as that defined by 150GeV < Ami < 200GeV and 95GeV < Am 2 < 120GeV for all panels in 
Figj2j Consequently, the boundaries of the CMS experiment in the cases of rriq = 650, 700GeV are mainly 
determined by the SR-II, while in the case of m, > 750GeV, they are determined by the SR-III. 
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• For squarks lighter than about 1200GeV, the tightest constraint comes from the CMS 
experiment, which is similar to the heavy gluino case; but with the increase of the 
squark mass, the strongest constraint may come from the ATLAS preliminary search 
for jets + Etp iss signal, which is shown in the lower right panel of Fig j3} The underlying 
reason for such a feature is that the capabilities of the SRs defined in the experiments 
in limiting SUSY depend on the configuration of the SUSY spectrum, and without 
the specification of the spectrum, there is no definite conclusion about which is the 
strongest. To illustrate this, we take one point from each of the two lower panels in 
Fig|3] as an example, and show the details of the constraints in Table [T| This table 
indicates that for the points Si and S 2 , the tightest constraints of the ATLAS prelim¬ 
inary search for jets + E™ lss are the SRs CM and BM respectively, and those of the 
ATLAS updated search correspond to the SRs 2jm and 2jt respectively. For the point 

51, the constraint from the former experimental analysis is weaker, while for the point 

5 2 , the situation reverses. 

• We checked that the ATLAS excess can be explained at 2cr level by squarks with mass 
up to about 1.4TeV. The wider mass range in comparison with the heavy gluino case is 
mainly due to the larger rate of the squark pair production in light gluino case. We also 
checked that the central value of the excess can be achieved for m g - < 1.2TeV if only 
the constraints from the ATLAS searches for the jets + E™ lss signal are considered, 
and that if the constraint from the CMS experiment is further considered, the maximal 
event number is only 9.5. The latter fact implies that the heavy gluino case is able to 
provide a slightly better explanation. 

• Note that for the lower right panel where rrig = 1.36TeV and nig = 1.5TeV, the effect 
from the squark-gluino associated production on the LHC searches is non-negligible. 
Discussing such an effect is beyond the scope of this work. 

Before we end this section, we have four comments about our explanation. First, com¬ 
paring Fig|2]and Figj3] in this work with Fig.2 in [5j where gluino pair production with the 
decay mode g —> qq \° —* qqZ\\ was used to explain the excess, we conclude that the squark 
explanation allows for a significantly wider range of m g - as well as a less compressed SUSY 
mass spectrum. One underlying reason we think is that the cut efficiency of the ATLAS 
on -Z search is usually larger for the squark explanation than for the gluino explanation, and 
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consequently, moderately heavy squarks are still capable of explaining the excess. Second, 
we emphasize again that in our simple scenario, both the event number of the ATLAS signal 
and the constraints are determined by few SUSY parameters. Consequently, the capability 
of our scenario to interpret the excess is limited. As we mentioned at the end of Sec.11, a 
more complex scenario with higher dimensional SUSY parameters may improve this situ¬ 
ation, but it needs tremendous calculation to search for the relevant parameter space pUj- 
Third, we note that so far there have appeared other simple SUSY scenarios to explain the 
excess, and the studies of these scenarios indicated that they can explain the excess at la 
level in certain narrow SUSY parameter spaces mm- This conclusion is slightly better 
than ours, where the best explanation is about 1.2a away from the central value of the 
excess. Three factors may contribute to the difference: 

• The difference in the theoretical hypothesis on SUSY, which determines the kinemat- 
ical distributions of the SUSY signals. For example, both HU and ra utilized the 
production pp — > gg with the loop-induced decay g — > gy§ — > gZ\ (x denotes the 
lightest neutralino in HU and the gra.vitino in [12]) to explain the excess. Comparing 
their interpretations with ours, one can learn that, although all of them considered 
the two-body decay of a strongly produced SUSY particle, due to the difference of 
the properties for the parent sparticles such as their spins and production channels, 
their kinematical distributions may differ greatly even when their predictions on the 
event number of the excess are same. One can get this point by comparing the Ej' mss 
distribution of the benchmark points P\ and P 2 in this paper with that of the best 
point in [12], which are presented in Figjijof this work and in Fig.6 of [2] respectively. 
As a result of the difference, in principle there might exist SUSY points for which the 
ATLAS leptonic-Z signal is moderately enhanced and meanwhile the CMS signal is 
appropriately suppressed. 

• The uncertainties induced by related simulations. For all scenarios to explain the 
excess, simulations of the experimental searches for SUSY have to be done. As shown 
in the appendices of [8] and where the validations of the simulations were explicitly 
presented, the uncertainties of the simulations are at the level of 20%, and therefore, 
the calculation performed by different groups may result in a significant deviation. As 
far as our simulations are concerned, the computed efficiency for the ATLAS signal 
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event is less than that presented by the ATLAS collaboration for the selected SUSY 
point by about 10%, and by contrast our efficiency for the CMS search is slightly larger 
than that in the CMS report. 

• The treatment of the CMS constraint. From the CMS report presented in [20j, one 
can only infer the approximate value of the S^ s for the SR-III. Confronted with such 
a situation, we in our previous work [8] calculated the S ^ s s for all the six SRs by the 
asymptotic CL S prescription |28j (see Table 2 of (8]). Furthermore, we pointed out 
that the SR-II may provide a stronger constraint on the parameters than the other 
SRs in discussing Fig|2] of this work. This conclusion indicates that the calculation of 
all Sj £ s s is necessary; but on the hand, since the values of the S^ s s were not explicitly 
given in other previous literatures and meanwhile they depend on calculation method, 
there might exist deviations for different authors in considering the CMS constraint. 
Let’s show this point by an explicit example. During the revision of this manuscript, 
the paper m appeared to interpret the excess in the NMSSM extension with a Dirac 
gluino, and the authors presented the details about their calculation of S^ s s. Briefly 
speaking, the calculation in [29j differs from ours in at least two aspects. One is the 
authors of [29] used the standard Bayesian procedure in the calculation, while we used 
the asymptotic CL S method [28]. The other is that the work [29] had considered the 
theoretical uncertainty in calculating the signal, and by contrast, we ignored such an 
effect. As a result of the differences, the values of S^ s s in ra are usually larger than 
our predictions by about 15%, and consequently, the CMS constraint is significantly 
relaxed in [29 ]. 

Finally, we note that in our scenario the singlino-dominated LSP is usually heavier than 
about 450GeV, and one may wonder how such a heavy dark matter (DM) achieves its 
measured relic density. In this case, the possible annihilation final states of the DM include 
//, VV, HiHj, AiAj and HiAj, where / (V) denotes any of the fermions (vector bosons) 
in the SM, and Hi ( Aj ) denotes a CP-even (CP-odd) Higgs boson (see [H] and references 
therein). The easiest way to achieve the density is through the s-channel annihilations 
mediated by a singlet-dominated Higgs boson, where just like the light DM case discussed 
in P2], the Higgs boson mass as well as the self-coupling coefficient k for the singlet fields 
play an important role in tuning the annihilation rate. 
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Pi 

P 2 

P 3 

SM 

7 rig (GeV) 

4500 

4500 

750 

- 

rriq (GeV) 

650 

700 

4500 

- 

m^o (GeV) 

530 

500 

650 

- 

m^o (GeV) 

430 

370 

565 

- 

Events in ATLAS on-Z SR(8TeV,20.3/6~ 1 ) 

11.0 

18.2 

10.5 

10.6 

X 2 E mis3 (d.O.f. = 9) 

7.6 

5.6 

8.0 

14.8 

XH T ( d -°-f- = 7 ) 

5.7 

2.6 

6.0 

13.8 

Xrtj {d.o.f. = 5) 

6.9 

4.7 

6.8 

14.1 

Xtotal (d-O.f. = 21) 

20.2 

12.9 

20.7 

42.7 

Events in ATLAS on-Z SR(14TeV,10/&- 1 ) 

172.8 

168.2 

103.6 

18.2 

Expected significance (14TeV, 10 fb ~ 1 ) 

25 

24 

14 

- 

Expected significance (8TeV, 20.3/6~ 1 ) 

2.4 

4.0 

2.3 

- 


TABLE II: Detailed information of the benchmark points Pi, P 2 and P 3 , including their contribu¬ 
tions to the ATLAS excess and behaviors at the LHC-14. 

IV. DISTRIBUTIONS OF THE EXCESS 

In this section, we investigate whether our explanation can reproduce the distributions 
of the excess reported by the ATLAS collaboration. For this end, we concentrate on three 
benchmark points PI, P2 and P3 with point P3 taken from our previous work [8]. Point 
PI and P3 correspond to the best points after considering all the constraints in the squark 
explanation and the gluino explanation respectively, and contribute 11 and 10.5 events to 
the excess. By contrast, point P2 only satisfies the constraints from the ATLAS jets +E™ lss 
searches, but it can reproduce the central value of the excess. Detailed information of these 
points is presented in Table [TTJ 

In order to compare our explanation with the experimental data for various distributions, 
we generate the distributions of E™ lss , Ht (the scalar sum of the Pts for the leptons and 
signal jets) and the jet multiplicity rij in the electron and muon combined channel for each 
parameter point. I 11 getting the distributions of E™ lss and H E , we include the overflow 
events into the last bin. The corresponding results are shown in Figj4| where the black solid 
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Ey W5 [GeV] i?T[GeV] Jet Multiplicity 


FIG. 4: Distributions of E™ lss , Ht and the jet multiplicity predicted by the three benchmark 
points. As a comparison, the experimental data for the ATLAS on-Z excess with the expected SM 
background extracted are also shown. Note that as did in [5], only the statistical uncertainties are 
included in the error bars of the data. 


circle with the error bar stands for the data obtained by the ATLAS experiment with the 
expected SM background subtracted [2], and the predictions of the points PI, P2 and P3 are 
marked by triangles, squares and asterisks respectively. To quantize the difference between 
the theoretical predictions and the corresponding experimental data, we define a y 2 function 
for each distribution by a simple way 


X 


2 


y ( Si - ^ ; ) 2 

V (^) 2 


(3) 


where .s* is the theoretical prediction in the ith bin, s t is the corresponding experimental 
datum, and 8s r is the error of the datum. In Table pd) we show the values of y 2 for different 
distributions. 

From Figj4]and Table |ITj one can learn that all points, especially point P2, can reproduce 
the distributions in an excellent way. The only significant difference between the squark 
explanation and the gluino explanation in generating the distributions comes from the jet 
multiplicity, i.e. for the former explanation, the rij distribution peaks at 3, while for the 
latter explanation, it peaks at 4. Due to the large errors of the data at current time, we can 
not determine which explanation is preferred to account for the excess. 
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V. TEST OF OUR EXPLANATION AT THE 14 TEV LHC 


Considering that the squark pair production rate at the LHC-14 can be greatly enhanced 
in comparison with that at the LHC-8, one may expect that the squark explanation will be 
tested very soon at the LHC. We investigate this issue by considering the lepton-Z + jets-\- 
Ejp iss signal of the production at the LHC-14. For simplicity, we assume the same cuts as 
those of the ATLAS on -Z search at the LHC-8, and estimate the SM background of the 
signal. In doing this, we suppose that the dominant background at the LHC-14 comes from 
the same processes as those at the LHC-8, which include flavor-symmetric backgrounds, 
Z + jets, rare top and diboson|[2]. Since it is hard to get accurate background events by 
directly simulating the processes at the LHC-14, we simulate each background process at 
the LHC-14 and LHC-8 separately to get the ratio of their rates after cuts, then we scale 
the background at the LHC-8, which was given in the ATLAS report [2], by this ratio. We 
realize that the results obtained in this way may deviate significantly from their true values, 
but without any detailed information of the ATLAS detector at the LHC-14, our results 
may serve as a rough estimate of the background. 

Once we know the signal and the total background after the cuts, we can calculate the 
expected significance by following formulae 


S 


N s 

y/Nh + eW*’ 


(4) 


where N s and Ay denote the event numbers for the signal and the background respectively, 
and the coefficient e parameterizes the effect induced by systematic errors. In our calculation, 
we set e = 30% for 14TeV LHC, which was adopted at the LHC-8|2]. 

Assuming 10/6 -1 integrated luminosity at the LHC-14, we present in Table |TT] the predic¬ 
tions of the points on the event number of the signal as well as that of the background. This 
table indicates that the signals at the LHC-14 after cuts are enhanced by more than 9 times, 
and by contrast the background is enhance by only about 2 times. As a result, either the 
squarks/gluino predicted by the points will be discovered, or in case of the non-observation 
of the leptonic-Z + jets + E™ lss signal, the points will be excluded. Moreover, comparing 
the squark explanation with the gluino explanation, we note that the former is more readily 
tested at 14TeV LHC. 
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VI. CONCLUSION 


In this paper, we intended to explain the 3cr excess recently reported by the ATLAS 
collaboration in the search for the leptonic -Z+jets+E™ lss signal. For this end, we considered 
the pair production of the first two generation squarks in the NMSSM with the decay 
chain q —> qx® QXi^- I 11 order to maximize the signal rate and also to simplify our 
analysis, we considered a singlino-dominated Xi and a bino-dominated x®i an d assumed 
both Br(q — > q\ °) and Br(x 2 —> X°\Z) at roughly 100%. With these assumptions, the 
parameters relevant to our analysis include the common squark mass rriq, the gluino mass 
rrig as well as the mass splittings Ami = rriq — and Am 2 = m^o — 0 . 

In order to find the parameter space that can explain the excess, we fixed rrig = 4.5TeV 
and rrig = 1.5TeV separately, and for each case, we performed a grid scan over the parameters 
Ami and Am 2 by choosing a series of rriq. Then for each parameter point we encountered, 
we simulated the process pp —> qq —>• x\ZqXiZq with the cuts adopted by the ATLAS on-Z 
search, the CMS on-Z search, and the ATLAS 2 ~ 6 jets -\-E™ lss searches respectively. Based 
on our simulations, we have following conclusions: 

• After considering the constraints from the ATLAS searches for jets + Etp iss signals, 
the central value of the ATLAS Z-peaked excess can be obtained for rriq < 1.2TeV. 

• If the constraint from the CMS on-Z search is further considered, more than 10 signal 
events are still attainable for rriq < 750GeV. 

• For the squarks as heavy as about 1.4TeV, the squark pair production can still account 
for the excess at 2a level without conflicting with any constraints. 

• Compared with the explanation by gluino pair production, the squark explanation 
allows for a significantly wider range of rriq as well as a less compressed SLISY mass 
spectrum. 

Moreover, we also investigated whether the squark pair production can reproduce the 
distributions of the excess reported by the ATLAS collaboration. We found that, quite 
similar to the gluino pair production, the squark pair production can fit the data quite well. 
Probing the squark explanation at the 14 TeV LHC was also investigated, and we concluded 
that only with 10/6” 1 integrated luminosity, the squarks that are able to explain the excess 
will be either discovered or excluded. 
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